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We study the coupled even number of microcavities with the balanced gain and loss between any 
pair of their neighboring components. The effective non-Hermitian Hamiltonian for such structure 
has the cyclic permutation-time symmetry with respect to the cavity modes, and this symmetry 
determines the patterns of the dynamical evolutions of the cavity modes. The systems also have 
multiple exceptional points for the degeneracy of the existing supermodes, exhibiting the “phase 
transition” of system dynamics across these exceptional points. We illustrate the quantum dynamical 
properties of the systems with the evolutions of cavity photon numbers and correlation functions. 

Moreover, we demonstrate the effects of the quantum noises accompanying the amplification and 
dissipation of the cavity modes. The reciprocal light transportation predicted with the effective 
non-Hermitian models for the similar couplers is violated by the unavoidable quantum noises. 


I. INTRODUCTION 


Coupled photonic structures exhibit rich physics and 
have wide applications. One category that has attracted 
wide research interest in the past decade is with prop¬ 
erly balanced gain-loss profiles for the coupling optical 
modes, as it is regarded as a realization of optical ana¬ 
logue for parity-time (PT)-symmetric quantum mechan¬ 
ics 0 , [2|. These proposed systems include optical waveg¬ 
uides [3 jp] , optical lattices |l-@], microcavities [3, H] , and 
others 9|, lL0|. So far most of the researches on these 
systems are concerned with the properties such as light 
transportation, and some of the interesting features of 
the systems have been demonstrated by the experiments 
with different setups EMI- An important character of 
a PT-symmetric system is the existence of the thresh¬ 
old known as exceptional point 0, G3 , across which 
the optical modes undergo a “phase transition” between 
periodic oscillation and exponential growth/decay. This 
feature can be applied to control the transmitted light 
intensity and enhance the added nonlinearity for gener¬ 
ating nonclassical light field states with PT systems 0 ] . 

In addition to the above mentioned properties, the 
quantum features of linear PT systems have received at¬ 
tention recently B00- The quantum dynamics of 
these systems should be well understood in view of the 
potential applications. On the one hand, the dynamics of 
these systems has been studied with PT-symmetric non- 
Hermitian Hamiltonians, which have real and complex 
eigenvalues, respectively, across the exceptional point 
where there is the unit ratio of amplification (dissipation) 
rate over coupling strength. The non-Hermitian system 
Hamiltonians lead to unique properties of quantum state 
transfer in coupled PT-symmetric arrays S0. On the 
other hand, as in any open quantum system, quantum 
noises coexisting with the gain and loss of the optical 
modes will inevitably affect the dynamics of the systems. 
The associated noises destroy the nonclassicality of quan¬ 
tum light sent into a linear PT-symmetric coupler 0. 


The external quantum noise driving fields, as well as the 
amplification and dissipation of light fields, can be de¬ 
scribed in terms of a stochastic Hermitian Hamiltonian 
poj . A clearer understanding of the possible differences 
between the two pictures is meaningful to further re¬ 
searches on the quantum aspect of similar coupled gain- 
loss systems. 

In this paper we propose a generalized structure from 
PT symmetry. The concerned linear quantum systems 
have alternately distributed equal gain and loss for their 
neighboring microcavity modes, so that their dynami¬ 
cal equations without external driving terms are invari¬ 
ant under cyclic permutation operations plus time re¬ 
versal operation for the coupling modes. One conse¬ 
quence of this cyclic permutation-time (CPT) symmetry 
is the existence of multiple exceptional points. We will 
illustrate the dynamical behaviors of the systems across 
these exceptional points. The other purpose of the cur¬ 
rent work is to understand the effects of quantum noises 
on this structure. The light transportation through a 
CPT system is compared for both non-Hermitian model 
and stochastic Hermitian Hamiltonian including quan¬ 
tum noises. We show that the corrections from the noises 
are considerable to such linear dynamical systems if the 
transmitted light is not so strong. 

The rest of the paper is organized as follows. In Sec. II 
we provide the description of the CPT systems with the 
non-Hermitian effective Hamiltonian, as well as with the 
stochastic Hermitian Hamiltonian including the quantum 
noises. The quantum dynamics of the systems is dis¬ 
cussed with an example in Sec. Ill; here we illustrate 
the evolution of cavity photon numbers and a kind of 
quantum correlation functions in the different regimes 
across the exceptional points. In Sec. IV we study the 
quantum noise effects to show their influence on the light 
transportation through a CPT system. The difference 
between the non-Hermitian Hamiltonian approach and 
the stochastic Hermitian Hamiltonian approach will be 
discussed in this section. Finally we conclude the cur- 
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FIG. 1: (color online) A schematic of CPT systems with even 
number of microcavities, which carry the balanced gain and 
loss (at the rate k) for each pair of the neighboring com¬ 
ponents. These cavities are coupled with each other by an 
adjustable strength J. In the non-Hermitian model for the 
system, the light transported from Fiber 1 to Fiber 2 has the 
same intensity as in the reverse process. Here we also con¬ 
sider the quantum noises as the random driving fields into the 
cavities. 


rent study in Sec. V. 


II. CPT SYSTEMS AS NON-HERMTIAN 
MODEL AND UNDER QUANTUM NOISES 

A generalized CPT system from the PT-symmetric 
ones is illustrated in Fig. |T| Here the microcavities carry 
the balanced gain and loss at the rate k for each pair of 
neighboring components, which couple with the strength 
J. Without loss of generality we let the external coherent 
driving field act on one of the gain cavities. This config¬ 
uration of even numbers of microcavities with the bal¬ 
anced gain and loss possesses interesting patterns for the 
involved dynamical supermodes illustrated below, which 
do not exist in coupled systems of odd numbers of cav¬ 
ities. We use the mode operator di (i = l,2,---,n) to 
denote those of the amplification cavities, while those of 
the dissipation ones are denoted as bj. They satisfy the 
commutation relation [c,,cj] = Sij for c* = a*,6i. The 
dynamical equation for the system reads (Ti = 1) 


d\ = iEe lAt , without loss of generality, where A is the 
detuning of the driving laser frequency from the intrin¬ 
sic cavity frequency and E the amplitude. Similar to 
the non-Hermitian Hamiltonians in §1, 0, the dynami¬ 
cal process in Eq. ([lj corresponds to the effective Hamil¬ 
tonian 

H e ff = ina\di + • • • + i«a^d n 
— inb\bi — ••• — inblfin 
+ J(aib\ + a\bi) + •••-)- J(aibl L + a\b n ) 

+ iE(a\e iAt — d\e~ iAt ). (2) 

Without the external driving field, this Hamiltonian 
is invariant under the cyclic permutation operation 
CV(ai, a n , b n ) = (6i, 02 , • • •, b n , di) plus the time 

reversal operation Ti = —i, hence the symmetry of the 
dynamical equation Eq. ©. One consequence of such 
symmetry is the reciprocal light transportation through 
the systems, if added with an external driving field. This 
inherits the same property of a linear PT-symmetric cou¬ 
pler, i.e. the light driven into the active (gain) cavity and 
going out of the passive (loss) one has the same intensity 
as in the reverse process (into the passive one and out of 
the active one). 

In reality the amplification and dissipation of the cav¬ 
ity modes are also accompanied by the quantum noises. 
These noises act as the additional random driving fields 
to the coherent one in Eq. ©, and then the to¬ 
tal drive vector in the equation takes the form d(t) = 
(' iEe lAt , 0, • • •, 0) T + i\/2K(H (*), r)i(t), • ■ •,&(*), Vn{t)) T ■ 
The amplification noise operator £j(i) and the dissipa¬ 
tion noise operator rji(t) satisfy the correlation relation 
(ii(t)H(t')) = (r/i(t)r)J(t')) = 6(t-t'). The noises to¬ 
gether with the cavity gains will make contribution to 
the output photon numbers, known as spontaneous pho¬ 
ton generation [l9| . The gain and loss of the cavity modes 
in Eq. JT]) can be described by their coupling with the 
associated reservoirs, so the process is due to a stochas¬ 
tic Hamiltonian [20 ] (the notation for the amplification 
noise operator follows that in 0 ) 



m 


[M]c(t) + d(t), 

/ in J J \ 

J —IK 


J IK J 


V j 


J —IK J 


H — J {cLib\ + a\bi ) + •••-}- J (cl ib^ + a\b n ) 

+ iV2K{(a\£\(t) - ai£i {t)) + -b (aUnW ~ a n £ n (t))} 

+ iE(a[e iAt — a\e~ iAt ) 

+ iV2K{{b\f)i{t) - hfil(t)) + -b (bivnit) - b n f)l(t))}. 

( 3 ) 


in the interaction picture with respect to the self oscil¬ 
lation Hamiltonian of the microcavities, where c{t ) = 
(ai(t), bi(t), ■ ■ ■, a n (t), b n (t)) T , and d(t) denotes the ex¬ 
ternal drives. The external coherent driving filed is 
into the first gain cavity of such linear systems, i.e. 


Obviously the CPT symmetry will be broken for this 
stochastic Hamiltonian (excluding the coherent drive 
term). The consequence of the symmetry breaking for 
the linear couplers in Fig. 1 is one of the main points we 
will discuss below. 
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(A, m) 

N = 2 

(-Vac 2 -J 2 , l), (Vac 2 -J 2 , l) 

N = 4 

(—k, 1),(k, 1 ), (—Vac 2 - 4J 2 , 1 ), (Vac 2 — 4J 2 , 1 ) 

N = 6 

(-Vac 2 - J2, 2), - J 2 , 2), (-V« 2 - 4J2, 1), 


(Vk 2 -4J 2 ,1) 

00 

(-«, 1),(k, 1), (-Wc 2 - 2d 2 , 2), (V« 2 - 2 J 2 , 2), 


(-V« 2 - 4J2,1), (Vac 2 - 4J2, 1) 

N = 10 

(~\A 2 - (^)J 2 ,2), (^ 2 -(^I)J 2 ,2), 


( v / K 2(3^1) J2i2 ) ) (v / k 2 ( 3^ )J2]2)i 


(-Vac 2 - 4J2,1), (Vac 2 - 4J 2 ,1) 

N = 12 

(-K, !),(«, 1), (-Vac 2 - J2, 2), (Vac 2 - J\ 2), 


(-Vac 2 - 3J 2 , 2), (Vac 2 - 3J 2 , 2), 


(-Vac 2 - 4J2,1), (Vac 2 - 4J 2 ,1) 


TABLE I: Supermode spectra for CPT systems. Here we list 
the eigenvalues and their degeneracy of the matrices — i[M] for 
those with small cavity number N. A denotes the eigenvalues, 
and m is the degeneracy of the corresponding eigenvalue. 

III. DYNAMICS ACROSS MULTIPLE 
EXCEPTIONAL POINTS 

A. existence of multiple exceptional points 


The dynamics of the CPT systems is completely de¬ 
termined by the matrix [M] in Eq. ([1]). Unlike a Her- 
mitian matrix this CPT-symmetric matrix does not have 
a real spectrum decomposition in terms of the orthonor¬ 
mal eigenvectors as the supermodes. However, various 
interesting properties arise from its non-Hermitianity. 

The spectrum and supermodes of the matrix [M] have 
regular patterns with the growth of the cavity number 
N, though they do not have the exact formulas. We 
list their proper forms in Tab. I for some CPT systems 
with small cavity number N. The eigenvalues of — i[M\ 
appear in pairs of opposite signs. Their distribution pat¬ 
tern is perceivable; the single eigenvalue pair with the 
values ±k appears periodically with the number N, and 
all other ones except for this pair and another one have 
the double degeneracy which affects the forms of the cor¬ 
responding supermodes (eigenvectors). The single eigen¬ 
value pair ±\/ ac 2 — J 2 for N = 2 is the one that exists in 
all linear PT-symmetric systems, and this pair of eigen¬ 
values degenerate at k = J as the exceptional point. For 
the higher numbers N, there are more such points at 
which the pairs of eigenvalues will degenerate further into 
the single ones. We call them the exceptional points for 
the CPT systems. For instance, two exceptional points, 
k = J and k = 2 J, exist for the CPT system of N = 6 . 
Across an exceptional point by adjusting the ratio J/ac, 
one pair of the eigenvalues takes the transition between 
the real and imaginary values, thus giving rise to different 
dynamical behaviors of the system. The dynamical evo¬ 
lution of the CPT systems can be therefore very different 
from the unitary evolution of a closed quantum system. 


B. evolution of observables 

The cavity modes evolve according to Eq. 0 to the 
following: 


Ci{t) 



(4) 


where [• • -]ij denote the matrix elements, and c( 0 ) = 
(Si, bi, ■ ■ ■ , a n , b n ) T . The interesting features of this lin¬ 
ear dynamical evolution come from the CPT symmetry 
of the dynamical matrix [M], The CPT-symmetric ma¬ 
trix can be diagonalized as —i[M] = [A'][A][If ] -1 (the 
diagonal elements of [A] are the eigenvalues in Tab. I), 
implying e _ d M b = [Aje^^A'] -1 . Due to the double de¬ 
generacy of the eigenvalues shown in Tab. 1, the expec¬ 
tation values of at least two amplification (dissipation) 
cavity modes will evolve in the same way for any ratio 
J/ac, given an external driving field into any of the cou¬ 
pled cavities. Take the system of N = 6 for example, the 
averages (a 2 {t)) and ( 03 (f)), as well as ( 61 (f)) and ( 63 (f)), 
are the same under a coherent driving field into the first 
amplification cavity. 

The dynamics across the different exceptional points 
should be understood for the CPT systems. Two types 
of information, the control on the transmitted light in¬ 
tensity with the system parameters and the possible cor¬ 
relations between the output modes from the individual 
microcavities, are useful to their potential applications. 
The former is straightforwardly quantified by the cavity 
photon evolving inside the microcavities and, for the lat¬ 
ter, we are concerned with the phase coherence between 
the cavity modes c, here. The phase coherence, as an 
important quantum correlation, can be measured by the 
average contrast [ 2 ll - [ 2 Sj 


C Ci C 




2 1 (4 A?) 1 

(4 Ci) + ( c]cj) 


(5) 


in interference experiments, where <7 = di or bi. From 
the input-output relation, the measured output modes 
di in experiment are related to the cavity modes as di = 
y/Vrjci, where 7 is the coupling rate between the cavities 
and the optical fibers for outputting the light. There 
is an extra term for the cavity into which the external 
driving field is added. 

Some of the CPT systems with one pair of supermodes 
of the eigenvalue A = ±k (see Tab. I) are not well tun¬ 
able, because the exponentially growing one of them (in¬ 
dependent of the system parameter J ) will quickly dom¬ 
inate the dynamical evolutions of the observables. Here 
we demonstrate the evolution patterns of the cavity pho¬ 
ton numbers and average contrasts between the cavity 
modes with the example of N = 6 , which has two excep¬ 
tional points at k = J and k = 2 J, respectively. We as- 
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FIG. 2: (color online) Evolution of photon numbers and aver¬ 
age contrast in the regime k < J. Here the plots are obtained 
with J/k = 2.5, E/k = 20 and A = 0. The initial quantum 
state of the system is assumed to be a vacuum state, (a) illus¬ 
trates the evolution of (ojai) (the thinner purple curve) and 
(b\bi) (the thicker orange curve); (b) is about the evolution 
of the average contrast C a ^ between the two modes, (c) and 
(d) illustrate the corresponding quantities for the modes <22 
and 1 ) 2 . 


sume that the amplification/dissipation rate k in the sys¬ 
tem is fixed, and the coupling strength J between the cav¬ 
ities can be adjusted. We start from the strong coupling 
regime J > n. Then all eigenvalues for the matrix —i[M] 
are imaginary. In this regime the eigenvalues of the effec¬ 
tive non-Hermitian Hamiltonian in Eq. © are all real, 
corresponding to the commonly known PT-symmetric 
regime for N = 2. The cavity modes are simply oscil¬ 
lating there. Consequently the two quantities mentioned 
above evolve as the quasi periodic functions shown in Fig. 
[2j Compared with the PT-symmetric systems (the spe¬ 
cial case of N = 2 in our generalization), more harmonic 
modes contribute to these quasi periodic evolutions for 
a high cavity number N. For brevity we only pick out 
two pairs of cavity modes to illustrate the evolution of 
the concerned quantities with time, and those for the 
other modes are similar. The extreme situation in this 
regime is the strong coupling limit J ft, where the sys¬ 
tem behaves like a beam-splitter array used for photonic 
quantum information processing (see, e.g. l2lM2Gl] ). In 
this regime with J > n, the pure transportation of light 
from the externally driven cavity can dominantly deter¬ 
mine the photon numbers in other cavities, and the field 
intensities in them are the results of the effective inter¬ 
ferences from two identical “sources” (in the driven cav¬ 
ity) respectively transmitting light along the clockwise 
and counterclockwise path in Fig. 1, which are formed 
by a beam-splitter type coupling between the cavities. 
The diagonal cavity to the driven one is under a perfect 
constructive interference because the photons are trans¬ 
ported along two symmetric paths to it. In our example 
of driving the first amplification cavity, the second dissi- 






FIG. 3: (color online) Evolution of photon numbers and av¬ 
erage contrast in the regime J < k < 2J. In (a) the thinner 
(purple) solid curve stands for the evolution of (fcj&i) with 
J/k = 0.6 and the thicker (orange) one stands for that of 
(b\bi) with ,// k = 0.7. In (b) the curves stand for the evolu¬ 
tion of the average contrast between the modes fii and fei for 
the ratios in (a), respectively. Correspondingly the plots in 
(c) and (d) are about the evolution of ( a\a 2 ) [with the same 
J/k ratios as in (a)] and the average contrast between the 
modes a 2 and 62 , respectively. All other parameters are the 
same as in Fig. [ 2 ] 

pation cavity as its diagonal unit therefore has a higher 
intensity than its neighboring amplification ones [see Fig. 
He)], as the photons are transported along two paths on 
which there is a symmetry with = (b\b 3 (t)) and 

(a' 2 a 2 (t)) = (a\a 3 (t)). This phenomenon is actually a 
consequence of the CPT symmetry. 

Next we lower the coupling strength such that the sys¬ 
tem enters the regime J < k < 2 J between the two excep¬ 
tional points. The evolution of the cavity photon num- 



FIG. 4: (color online) Evolution of photon numbers and 
average contrast in the regime k > 2 J. All plots here describe 
the same quantities as in Fig. [3j except that the thinner 
(purple) curves now stands for the situation with J/k = 0.4 
and the thicker (orange) ones for J/k = 0.2. 
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bers and the average contrast between the cavity modes 
are shown in Fig. [3] In contrast to the situations in 
Fig. [21 the photon numbers significantly grow up with 
time, given the same external driving field as in Fig. [2] 
It is due to the transition of t wo oscilla ting supermodes 
(of the eigenvalue value A = y/n 2 — J 2 with the double 
degeneracy) to the exponentially increasing ones. These 
two supermodes will dominate in the end. As the con¬ 
sequence, the average contrasts in Fig. [3] will become 
stable with time. Comparing the results in Fig. [3ja) and 
[3] (c), one sees that the photon numbers in the amplifi¬ 
cation cavities begin to surpass those in the dissipation 
ones. As the coupling strength J is reduced, the pho¬ 
ton numbers in all cavities will steadily increase instead. 
The change of the average contrast with the coupling J 
shows the opposite tendency, due to the different growth 
paces for the correlation |(cjcj}| (i yf j) and the photon 
numbers. 

After the coupling strength J is lowered further to tune 
the system into the regime of k > 2 J, the evolutions of 
the concerned quantities illustrated in Fig. 0] still behave 
similarly to those in Fig. [3] but show the significant dif¬ 
ference in their quantitative values (in Fig. [4] we only 
show the beginning stage for the photon numbers to sig¬ 
nificantly deviate from those in Fig. [3]). A main character 
in this regime is the saturation of photon amplification 




FIG. 5: (color online) Comparison between the different dy¬ 
namical regimes for a CPT system of N = 6. (a) shows the 
photon number (6[&i) distributing with the parameter J/k 
at the time nt = 20, when the supermodes of the largest 
eigenvalue have dominated the system evolution, (b) is the 
corresponding distributions of the average contrast C ai b 1 at 
two different moments, under a driving field of E/k = 20. 
Here we mark the two boundaries of the three regimes with 
the dashed vertical lines. The photon numbers in other cavi¬ 
ties and the average contrasts between other modes have the 
similar distributions. 


with the decrease of the coupling intensity J ; see Fig. 
[3(a). The photon number growth rate determined by 
the factors e VK 2 -J 2 d--r) and 6 Vk 2 -4 J 2 (t-r) f rom ex _ 
ponentially increasing modes gets larger for a smaller J. 
However, the smaller J will reduce the amount of pho¬ 
tons transported from the driven cavity to the rest ones. 
These two competing tendencies result in the maximum 
cavity photon numbers at a coupling strength J inside the 
regime. In dynamics this regime differs from the other 
one illu strated in Fig. [3] by having one more supermode 
(A = y/ k 2 — 4J 2 ) to undergo a transition from periodi¬ 
cally oscillating to exponentially growing. The distinc¬ 
tion of the two regimes possessing growing supermodes 
can be better seen in Fig. Ob), which shows an aver¬ 
age contrast distribution after a short and longer period 
of time. Though the supermodes with the largest eigen¬ 
value A = y/ k 2 — J 2 will dominate nonetheless, the extra 
growing supermode in the regime k > 2J considerably 
affects the transient behavior of the correlation function, 
making it take longer time to become stable. In contrast, 
the correlation function as the average contrast quickly 
stabilizes in the regime J < k < 2 J, causing the two dis¬ 
tribution curves in Fig. [3(b) mostly overlapped in this 
middle regime. It should be noted that, at any moment, 
all observables distribute continuously across the excep¬ 
tional points as seen from Fig. [3 This is totally different 
from a phase transition of many-body system, which is 
characterized by discontinuous change of system entropy 
or its derivatives in parameter space. Instead the dy¬ 
namical behavior of a PT or CPT system takes gradual 
transitions across the exceptional points, and the nature 
of such transitions is the passage of the supermodes be¬ 
tween their imaginary and real eigenvalues. 

Like beam-splitter arrays as an important in¬ 

gredient for linear optical quantum information process¬ 
ing [23], our concerned CPT systems, also as a type of 
development from photonic molecule structures [2a, [ 23 ], 
could be suitable for the similar and other applications. 
Here we have demonstrated the evolution of two sets of 
important quantities, the cavity photon numbers and the 
correlation functions (average contrasts). The unnormal- 
ized correlation functions (c\cj) from the latter constitute 
the correlation matrix or covariance matrix for a CPT 
system. Together with the photon numbers and aver¬ 
age cavity modes (c,), these quantities fully depict the 
Gaussian states of such linear quantum systems [30, (3l|. 

In the above discussion we have assumed a constant 
amplification rate to neglect the gain saturation for not 
too strong cavity fields. When the cavity field intensi¬ 
ties are sufficiently large, the amplification rate will be 
reduced with the growing cavity field intensities /,; pro¬ 
portional to (ajdj), i.e. it becomes re/(l + Ii/I s ) where 
I s is the saturation intensity. This will explicitly break 
the CPT symmetry of the concerned systems due to the 
loss of the balanced amplification and dissipation for the 
neighboring units. The proper operation of the systems 
should be under the condition Ii/I s <C 1. This condi¬ 
tion of weak cavity fields compared with the saturation 
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FIG. 6 : (color online) Difference between transported photon numbers along the mutually reverse paths. The solid curves 
stand for the transported photon number differences, which are proportional to the difference in the measured light powers. 
Here we illustrate the quantity 1 ( 6363 )^ — (a{ai)s|, where ( 6363 )^ is the photon number of the third dissipation cavity, which is 
found in the forward transmission pattern with the coherent drive on the first amplification cavity, and (o}ai)B is the photon 
number in the first amplification cavity due to the backward transmission pattern of driving the third dissipation cavity. Such 
difference comes from the amplification noises. The dashed curves represent the constant forwardly and backwardly transported 
photon numbers purely due to a coherent driving field with the intensity E/k = 5, which are found with the non-Hermitian 
model, (a), (b) and (c) are about the ratio J/k = 1.2, 0.6 and 0.4, respectively. 


value necessitates the consideration of quantum noise ef¬ 
fect discussed below. 


IV. INFLUENCE OF QUANTUM NOISES ON 
LIGHT TRANSPORTATION 


Due to CPT symmetry, the dynamical evolution ma¬ 
trix in Eq. Q is a symmetric one. If the systems 

are only described in terms of the non-Hermitian model 
in Eq. ©, one consequence of the symmetry is recipro¬ 
cal light transportation. The light input by an external 
driving field into one amplification microcavity and out 
of another dissipation one would have the same intensity 
as in the reverse process (into the latter and out of the 
former). The violation of this property, known as opti¬ 
cal isolation, from the saturation nonlinearity of the gain 
process in a PT system has been experimentally demon¬ 
strated recently [15]. Here we will illustrate a similar 
phenomenon of different nature. 

The amplification and dissipation of the cavity modes 
are inevitably accompanied by the associated quantum 
noises. Generally these quantum noises are neglected in 
all studies of the similar systems except in 18, 19]. The 
quantum noises in the CPT systems behave as the ex¬ 
ternal random driving fields in Eq. ©, in addition to 
the coherent one to implement the light transportation. 
In a PT system the noise accompanying the gain pro¬ 
cess makes extra contribution to the photon numbers of 
the coupling modes 0. For the CPT structure gener¬ 
alized from PT symmetry, there are more cavities under 
gain noise, and these noises act independently but af¬ 
fect the photon numbers in all coupled cavities. In an 
actual light transportation, therefore, the amplification 
quantum noises as the external random drives should 
also be considered in a strict sense. Moreover, the noise 
driving field intensities are independent of the coherent 
one, and are simply determined by the fixed amplifica¬ 
tion/dissipation rate n. Though two coherent driving 


fields with the same intensity transmit symmetrically 
along the mutually reverse paths between an amplifi¬ 
cation and a dissipation cavity, the quantum noises of 
amplification add asymmetric contributions to this situ¬ 
ation, leading to the violation of reciprocal light trans¬ 
portation from the non-Hermitian model. 

We still work with the system of N = 6 for il¬ 
lustration. The light transportation can be between 
the first amplification cavity and the third dissipation 
cavity as an example. The non-reciprocal light trans¬ 
portation under the quantum noises can be understood 
from the response of the system to the external driv¬ 
ing fields. In the non-Hermitian model, the recip¬ 
rocal transportation from a coherent driving field is 
due to the symmetric matrix elements = 

[ e -i[M](t-r)] 6 i gj The amplification noises ef¬ 

fect on the first gain cavity through the matrix ele¬ 
ments [e -i [ M K t- ' r )]n, [ e -ffM](t-r)]i 3 and [ e -i[M](i-r)] i5) 

while drive the third loss cavity through the elements 
t)] 61j [ e -i[M](t-r)] 63 and [ e -i[M](i—r)] 65 . Thege 

two sets of matrix elements have no symmetric relation. 

In Fig. [G] we illustrate the difference in the mutually 
reverse light transportations for the three regimes across 
the exceptional points. These differences in the cavity 
photon numbers are proportional to the difference of the 
measured light powers. Here we also compare the results 
with those obtained with the non-Hermitian model for a 
transmitted light of moderate intensity, showing that, in 
some of the regimes, the corrections from the quantum 
noises can be even more significant than the transmit¬ 
ted light itself. Similarly the quantum noises also affect 
the correlation functions we have discussed before. The 
quantum noise effects will become more prominent in the 
CPT systems of higher cavity number N. Since the back¬ 
ground noises act as the random driving fields with the 
intensity proportional to \J~k [see Eq. (J3J)] , their effects 
will be overshadowed by a coherent one of sufficiently 
high intensity E n (this value should go up with the 
cavity number N). However, the cavity field intensities 
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should be limited to avoid the significant gain saturation 
that breaks the balanced gain and loss between the neigh¬ 
boring units. The amplification noises could contribute 
to considerable amount of the transported photons in 
addition to those from a not too strong coherent driv¬ 
ing field, and should be taken into account in the general 
performance of the systems. 

The non-reciprocal light transportation is an appro¬ 
priate phenomenon to show the non-equivalence between 
the non-Hermitian model in Eq. © and the more realis¬ 
tic description of the similar systems with Eq. © . Con¬ 
sidering the quantum noise effects, the concerned systems 
do not exactly possess PT or generalized CPT symme¬ 
try. However, one of the the most important features for 
the systems is the dynamical transitions across the ex¬ 
ceptional points. The existence of the exceptional points 
can be regarded as the consequence of a “partial sym¬ 
metry” , i.e. the PT or CPT symmetry of the dynamical 
matrix [M] in equations of motion. 

V. CONCLUSION 

We have studied a generalized structure from the 
widely explored PT symmetric optical couplers. Due 


to the CPT symmetry of the dynamical matrix of such 
systems, they possess multiple thresholds (exceptional 
points) for the transitions of system dynamics. The dy¬ 
namics of an example system in the regimes separated 
by its exceptional points is illustrated with the evolution 
of cavity photon numbers and correlation functions. We 
find that the exceptional points of the CPT systems dif¬ 
fer from the transitional points of many-body systems by 
nature; no clear-cut boundary exists for the transitions of 
physical quantities across them. On the other hand, we 
demonstrate the importance of the quantum noises with 
their effect on light transportation. The quantum noise 
effects should be considered when the external coherent 
driving fields may not be very strong. The understand¬ 
ing of the dynamics of such systems is meaningful to their 
potential applications. 
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